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Abstract 

In  an  earlier  paper,  the  authors  derived  a Fredholm 
integral  equation  of  the  first  kind  for  the  solution  of  the  inverse 
source  problem  for  acoustic  waves.  The  eigenvalues  of  this  equation 
were  shown  to  converge  rapidly  to  zero  and  also  to  include  zero. 
Thus,  the  solution  was  shown  to  be  non-unique  and  even  the  parti- 
cular part  of  the  solution  of  that  equation  was  ill-conditioned.  In 
this  note  it  is  shown  how  to  obtain  the  non-trivial  information  of 
that  integral  equation  in  a well -conditioned  manner. 


In  an  oarlier  paper the  authors  developed  analytical 
characterizations  of  the  non-uniqueness  of  solutions  of  the  inverse 
source  problem  in  acoustics  and  electromagnetics.  In  this  problem 
the  field  radiated  by  a source  distribution  is  observed  on  a closed 
surface  (e.g.,  a sphere)  outside  the  region  containing  tne  source. 

The  objective  is  to  obtain  information  about  the  source. 

It  was  further  shown  in  Ref.  1,  that  the  source  distribution 
satisfies  a Fredholm  integral  equation  of  the  first  kind  with  eigen- 

Of\  ^ 

values  which  rapidly  approached  zero  0 (n  ) in  addition  to 
the  eigenvalue  zero,  itself.  Thus,  even  the  -gartiaulca?  solution  of 
this  equation— for  that  part  of  the  source  which  could  be  determined 
from  the  radiated  field — is  ill-oonditioned\  i.e. , highly  unstable 
to  noise  in  the  "higher"  eigenfunctions. 

The  purpose  of  the  note  is  to  show  that  the  information 
contained  in  that  integral  equation  can,  in  fact,  be  obtained  in  a 
well-conditioned  manner.  The  discussion  here  will  be  limited  to  the 
acoustic  case;  the  extension  to  the  electromagnetic  case  is 
straightforward. 

It  is  assumed  that  u(x,  w),  x = (x  , x , x ),  is  a 
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solution  of  the  inhomogeneous  Fourier  (time)  transformed  wave 
equation. 


(V^  + w^c”^)  u(x,  w)  = -f(x,  (d). 


(1) 


subject  to  the  radiation  condition 
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u(x,  toK  u (x,  w)  exp{iwx/c}  /(4tt  x),  x ->  » . (2) 

“ 0 

Here,  x denotes  a unit  vector  in  the  direction  of  and  x = |x|  is 
the  magnitude  of  the  vector  _x. 

The  source  function  f()^,  u)  is  assumed  to  be  confined  to 
the  interior  of  a sphere  of  radius  a.  For  observations  outside  of 
this  sphere,  the  solution  has  the  integral  representation 

u(x,  w)  = iwc"*  f)  E 
m=0  n=  -m 


Here, 

a 

0 


(4) 


and  are  respectively,  the  spherical  Bessel  function  and 
Hankel  function  of  the  first  kind  and  is  the  spherical  harmonic 
of  order  mn.  The  functions  ^n,j,{x,  w)  are  defined  by 


f™„(x,  w) 
mn 


IT  Zir 

= j siMde  / d(j)  f(x,  a))Y^*  (e,(})). 


(5) 


These  functions  are  the  coefficients  of  f in  its  spherical  harmonic 
expansion 


(6) 


“>  m 

f(>L.  “)  = Z Z (X. 

n=o  n=  -m 

The  representation  (3)  arises  from  the  Green's  function 
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representation  of  the  solution  and  the  spherical  harmonic 

expansion  of  the  Green's  function  in  Ref.  2 on  p.  742. 

Since  the  spherical  harmonics  are  a complete  set  of 

functions,  knowledge  of  the  coefficients  f^^^  constitutes  knowledge 

of  f itself.  However,  from  (3,  4,  5)  it  is  seen  that  the  radiated 

field  (i.e.,  u for  x > a)  is  a function,  not  o-''  but  of  their 

projeationa  c^j^^  on  the  spherical  Bessel  functioii  . Since  the 

set  of  functions  are  not  complete,  this  led  in  Ref.  i to 

one  analytical  characterization  of  the  non-uniqueness  of  the  solution 

of  the  inverse  source  problem. 

In  the  particular  solution  of  the  inverse  source  problem 

in  terms  of  spherical  harmonics,  one  could  only  hope  to  determine 

the  coefficients  c„„  and  not  the  coefficients  f„„.  It  will  now 
mn  mn 

be  shown  how  this  can  be  done  without  resorting  to  an  ill-conditioned 
integral  equation. 

To  begin,  the  function  \/[x,  u)  is  introduced,  denoting  any 
solution  of  the  homogeneous  reduced  wave  equation  in  the  region  0, 
on  the  boundary  of  which  (3D)  the  radiated  field,  u(x.,  w)  is  to  be 
observed.  Now  Green's  theorem  is  applied  on  the  domain  D to  the 
quantity 

v(2L,  a))[V^  + u)^c"^3u(x,  u)  - u(x,  a))[V^  + w^c”^]  v(x,  to). 

The  result  is 
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j v(x,  0,)  t(x,  »)  dV  = j {u^  - « . (7) 

X < a 3D 

For  any  function  v(x,  w)  of  the  prescribed  type,  this  is  an  integral 
equation  for  f(x,  w).  A particular  choice  suggested  by  the  discussion 
above  is  any  function  in  the  class 

'^mn  " m = 1,  2 |n|  < m . (8) 

With  v(x_,  w)  = 4'*  , (7)  becomes 


9D 


Thus,  the  are  given  in  terms  of  well -conditioned  operations— 
multiplication  by  known  functions  and  integration— on  the  observed 
data. 

In  Ref.  1,  an  equivalent  characterization  of  non-uniqueness 
was  given  in  terms  of  the  space-time  transform,  f{Jc,  w),  of  the 
source  distribution.  In  particular,  it  was  shown  that  the  radiated 
field  was  given  totally  in  terms  of  the  value  of  f(jk,  w)  on  the 
hyper-cone  where 
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k = lojJ/c  . 


(10) 


Conversely,  from  observations  of  the  radiated  field  and  its  normal 
derivative  one  can  only  determine  f{Jc,  w)  on  this  hyper-cone. 

These  function  values  can  also  be  extracted  from  (7)  in 
a well -conditioned  manner.  To  do  so,  one  need  only  choose  v to  be 
in  the  continuum  of  functions  over  unit  vectors  k. 


/V 

v()^,  0),  k)  = expCiwk  • x/c)  . 


(11) 


Now  (7)  becomes 


f(J(,  w) 


n u - exp{ijc  • dA, 


(12) 


A 

Here  n is  the  unit  outward  normal  to  9D. 
In  Reference  1,  the  function 


v(x,  w,  i)  = j^(w  R/c),  R = |x  - il  (13) 

was  used.  It  was  this  function  which  led  to  an  ill-conditioned 
integral  equation  for  f.  It  is  the  authors'  present  point  of  view 
that  this  is  to  be  rejected  as  a poor  choice  for  v(x,  w). 
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